A Fast Algorithm for Robust Regression with Penalised Trimmed Squares

The presence of groups containing high leverage outliers makes linear regression a difficult problem due to the masking effect. The available high breakdown estimators based on Least Trimmed Squares often do not succeed in detecting masked high leverage outliers in finite samples. An alternative to the LTS estimator, called Penalised Trimmed Squares (PTS) estimator, was introduced by the authors in \cite{ZiouAv:05,ZiAvPi:07} and it appears to be less sensitive to the masking problem. This estimator is defined by a Quadratic Mixed Integer Programming (QMIP) problem, where in the objective function a penalty cost for each observation is included which serves as an upper bound on the residual error for any feasible regression line. Since the PTS does not require presetting the number of outliers to delete from the data set, it has better efficiency with respect to other estimators. However, due to the high computational complexity of the resulting QMIP problem, exact solutions for moderately large regression problems is infeasible. In this paper we further establish the theoretical properties of the PTS estimator, such as high breakdown and efficiency, and propose an approximate algorithm called Fast-PTS to compute the PTS estimator for large data sets efficiently. Extensive computational experiments on sets of benchmark instances with varying degrees of outlier contamination, indicate that the proposed algorithm performs well in identifying groups of high leverage outliers in reasonable computational time.Comment: 27 page

Agent based modelling helps in understanding the rules by which fibroblasts support keratinocyte colony formation

Background: Autologous keratincoytes are routinely expanded using irradiated mouse fibroblasts and bovine serum for clinical use. With growing concerns about the safety of these xenobiotic materials, it is desirable to culture keratinocytes in media without animal derived products. An improved understanding of epithelial/mesenchymal interactions could assist in this. Methodology/Principal Findings: A keratincyte/fibroblast o-culture model was developed by extending an agent-based keratinocyte colony formation model to include the response of keratinocytes to both fibroblasts and serum. The model was validated by comparison of the in virtuo and in vitro multicellular behaviour of keratinocytes and fibroblasts in single and co-culture in Greens medium. To test the robustness of the model, several properties of the fibroblasts were changed to investigate their influence on the multicellular morphogenesis of keratinocyes and fibroblasts. The model was then used to generate hypotheses to explore the interactions of both proliferative and growth arrested fibroblasts with keratinocytes. The key predictions arising from the model which were confirmed by in vitro experiments were that 1) the ratio of fibroblasts to keratinocytes would critically influence keratinocyte colony expansion, 2) this ratio needed to be optimum at the beginning of the co-culture, 3) proliferative fibroblasts would be more effective than irradiated cells in expanding keratinocytes and 4) in the presence of an adequate number of fibroblasts, keratinocyte expansion would be independent of serum. Conclusions: A closely associated computational and biological approach is a powerful tool for understanding complex biological systems such as the interactions between keratinocytes and fibroblasts. The key outcome of this study is the finding that the early addition of a critical ratio of proliferative fibroblasts can give rapid keratinocyte expansion without the use of irradiated mouse fibroblasts and bovine serum

A Markov computer simulation model of the economics of neuromuscular blockade in patients with acute respiratory distress syndrome

BACKGROUND: Management of acute respiratory distress syndrome (ARDS) in the intensive care unit (ICU) is clinically challenging and costly. Neuromuscular blocking agents may facilitate mechanical ventilation and improve oxygenation, but may result in prolonged recovery of neuromuscular function and acute quadriplegic myopathy syndrome (AQMS). The goal of this study was to address a hypothetical question via computer modeling: Would a reduction in intubation time of 6 hours and/or a reduction in the incidence of AQMS from 25% to 21%, provide enough benefit to justify a drug with an additional expenditure of $267 (the difference in acquisition cost between a generic and brand name neuromuscular blocker)? METHODS: The base case was a 55 year-old man in the ICU with ARDS who receives neuromuscular blockade for 3.5 days. A Markov model was designed with hypothetical patients in 1 of 6 mutually exclusive health states: ICU-intubated, ICU-extubated, hospital ward, long-term care, home, or death, over a period of 6 months. The net monetary benefit was computed. RESULTS: Our computer simulation modeling predicted the mean cost for ARDS patients receiving standard care for 6 months to be$62,238 (5% – 95% percentiles $42,259 –$83,766), with an overall 6-month mortality of 39%. Assuming a ceiling ratio of $35,000, even if a drug (that cost$267 more) hypothetically reduced AQMS from 25% to 21% and decreased intubation time by 6 hours, the net monetary benefit would only equal $137. CONCLUSION: ARDS patients receiving a neuromuscular blocker have a high mortality, and unpredictable outcome, which results in large variability in costs per case. If a patient dies, there is no benefit to any drug that reduces ventilation time or AQMS incidence. A prospective, randomized pharmacoeconomic study of neuromuscular blockers in the ICU to asses AQMS or intubation times is impractical because of the highly variable clinical course of patients with ARDS

Decreased aortic growth and middle aortic syndrome in patients with neuroblastoma after radiation therapy

Background: Long-term CT follow-up studies are required in pediatric patients who have received intraoperative radiation therapy (IORT) and external beam radiation therapy (EBRT) to assess vascular toxicities and to determine the exact complication rate. Objective: To analyze with CT the effects of radiation therapy (RT) on the growth of the aorta in neuroblastoma patients.Materials and methods Abdominal CT scans of 31 patients with intraabdominal neuroblastoma (stage II–IV), treated with RT (20 IORT±EBRT, 11 EBRT alone), were analyzed retrospectively. The diameter of the abdominal aorta was measured before and after RT. These data were compared to normal and predicted normal aortic diameters of children, according to the model of Fitzgerald, Donaldson and Poznanski (aortic diameter in centimeters = 0.844 + 0.0599 × age in years), and to the diameters of a control group of children who had not undergone RT. Statistical analyses for the primary aims were performed using the chi-squared test, t-test, Mann-Whitney test, nonparametric Wilcoxon matched-pairs test and analysis of variance for repeated measures. Clinical files and imaging studies were evaluated for signs of late vascular complications of neuroblastoma patients who had received RT. Results: The mean diameter before and after RT and the growth of the aorta were significantly lower than expected in patients with neuroblastoma (P<0.05 for each) and when compared to the growth in a control group with normal and nonirradiated aortas. Among the patients who had received RT, there was no difference due to the type of RT. Seven patients from the IORT±EBRT group developed vascular complications, which included hypertension (five), middle aortic syndrome (two), death due to mesenteric ischemia (one) and critical aortic stenosis, which required aortic bypass surgery (two).Conclusion Patients with neuroblastoma who had received RT showed impaired growth of the abdominal aorta. Significant long-term vascular complications occurred in seven patients who received IORT±EBRT. Thus, CT evaluation of patients with neuroblastoma who receive RT should include not only reports of changes in tumor extension, but also documentation of perfusion, and the size and growth of the aorta and its branches over time

Comparative determination of HIV-1 co-receptor tropism by Enhanced Sensitivity Trofile, gp120 V3-loop RNA and DNA genotyping

BACKGROUND: Trofile is the prospectively validated HIV-1 tropism assay. Its use is limited by high costs, long turn-around time, and inability to test patients with very low or undetectable viremia. We aimed at assessing the efficiency of population genotypic assays based on gp120 V3-loop sequencing for the determination of tropism in plasma viral RNA and in whole-blood viral DNA. Contemporary and follow-up plasma and whole-blood samples from patients undergoing tropism testing via the enhanced sensitivity Trofile (ESTA) were collected. Clinical and clonal geno2pheno[coreceptor] (G2P) models at 10% and at optimised 5.7% false positive rate cutoff were evaluated using viral DNA and RNA samples, compared against each other and ESTA, using Cohen's kappa, phylogenetic analysis, and area under the receiver operating characteristic (AUROC). RESULTS: Both clinical and clonal G2P (with different false positive rates) showed good performances in predicting the ESTA outcome (for V3 RNA-based clinical G2P at 10% false positive rate AUROC = 0.83, sensitivity = 90%, specificity = 75%). The rate of agreement between DNA- and RNA-based clinical G2P was fair (kappa = 0.74, p < 0.0001), and DNA-based clinical G2P accurately predicted the plasma ESTA (AUROC = 0.86). Significant differences in the viral populations were detected when comparing inter/intra patient diversity of viral DNA with RNA sequences. CONCLUSIONS: Plasma HIV RNA or whole-blood HIV DNA V3-loop sequencing interpreted with clinical G2P is cheap and can be a good surrogate for ESTA. Although there may be differences among viral RNA and DNA populations in the same host, DNA-based G2P may be used as an indication of viral tropism in patients with undetectable plasma viremia

Plasmalogen enrichment in exosomes secreted by a nematode parasite versus those derived from its mouse host: implications for exosome stability and biology

Extracellular vesicles (EVs) mediate communication between cells and organisms across all 3 kingdoms of life. Several reports have demonstrated that EVs can transfer molecules between phylogenetically diverse species and can be used by parasites to alter the properties of the host environment. Whilst the concept of vesicle secretion and uptake is broad reaching, the molecular composition of these complexes is expected to be diverse based on the physiology and environmental niche of different organisms. Exosomes are one class of EVs originally defined based on their endocytic origin, as these derive from multivesicular bodies that then fuse with the plasma membrane releasing them into the extracellular environment. The term exosome has also been used to describe any small EVs recovered by high-speed ultracentrifugation, irrespective of origin since this is not always well characterized. Here, we use comparative global lipidomic analysis to examine the composition of EVs, which we term exosomes, that are secreted by the gastrointestinal nematode, Heligmosomoides polygyrus, in relation to exosomes secreted by cells of its murine host. Ultra-performance liquid chromatography – tandem mass spectrometry (UPLC-MS/MS) analysis reveals a 9- to 62-fold enrichment of plasmalogens, as well as other classes of ether glycerophospholipids, along with a relative lack of cholesterol and sphingomyelin (SM) in the nematode exosomes compared with those secreted by murine cells. Biophysical analyses of the membrane dynamics of these exosomes demonstrate increased rigidity in those from the nematode, and parallel studies with synthetic vesicles support a role of plasmalogens in stabilizing the membrane structure. These results suggest that nematodes can maintain exosome membrane structure and integrity through increased plasmalogens, compensating for diminished levels of other lipids, including cholesterol and SM. This work also illuminates the prevalence of plasmalogens in some EVs, which has not been widely reported and could have implications for the biochemical or immunomodulatory properties of EVs. Further comparative analyses such as those described here will shed light on diversity in the molecular properties of EVs that enable them to function in cross-species communication

From Sea to Sea: Canada's Three Oceans of Biodiversity

Evaluating and understanding biodiversity in marine ecosystems are both necessary and challenging for conservation. This paper compiles and summarizes current knowledge of the diversity of marine taxa in Canada's three oceans while recognizing that this compilation is incomplete and will change in the future. That Canada has the longest coastline in the world and incorporates distinctly different biogeographic provinces and ecoregions (e.g., temperate through ice-covered areas) constrains this analysis. The taxonomic groups presented here include microbes, phytoplankton, macroalgae, zooplankton, benthic infauna, fishes, and marine mammals. The minimum number of species or taxa compiled here is 15,988 for the three Canadian oceans. However, this number clearly underestimates in several ways the total number of taxa present. First, there are significant gaps in the published literature. Second, the diversity of many habitats has not been compiled for all taxonomic groups (e.g., intertidal rocky shores, deep sea), and data compilations are based on short-term, directed research programs or longer-term monitoring activities with limited spatial resolution. Third, the biodiversity of large organisms is well known, but this is not true of smaller organisms. Finally, the greatest constraint on this summary is the willingness and capacity of those who collected the data to make it available to those interested in biodiversity meta-analyses. Confirmation of identities and intercomparison of studies are also constrained by the disturbing rate of decline in the number of taxonomists and systematists specializing on marine taxa in Canada. This decline is mostly the result of retirements of current specialists and to a lack of training and employment opportunities for new ones. Considering the difficulties encountered in compiling an overview of biogeographic data and the diversity of species or taxa in Canada's three oceans, this synthesis is intended to serve as a biodiversity baseline for a new program on marine biodiversity, the Canadian Healthy Ocean Network. A major effort needs to be undertaken to establish a complete baseline of Canadian marine biodiversity of all taxonomic groups, especially if we are to understand and conserve this part of Canada's natural heritage

Determinants of recovery from post-COVID-19 dyspnoea: analysis of UK prospective cohorts of hospitalised COVID-19 patients and community-based controls

Background The risk factors for recovery from COVID-19 dyspnoea are poorly understood. We investigated determinants of recovery from dyspnoea in adults with COVID-19 and compared these to determinants of recovery from non-COVID-19 dyspnoea. Methods We used data from two prospective cohort studies: PHOSP-COVID (patients hospitalised between March 2020 and April 2021 with COVID-19) and COVIDENCE UK (community cohort studied over the same time period). PHOSP-COVID data were collected during hospitalisation and at 5-month and 1-year follow-up visits. COVIDENCE UK data were obtained through baseline and monthly online questionnaires. Dyspnoea was measured in both cohorts with the Medical Research Council Dyspnoea Scale. We used multivariable logistic regression to identify determinants associated with a reduction in dyspnoea between 5-month and 1-year follow-up. Findings We included 990 PHOSP-COVID and 3309 COVIDENCE UK participants. We observed higher odds of improvement between 5-month and 1-year follow-up among PHOSP-COVID participants who were younger (odds ratio 1.02 per year, 95% CI 1.01–1.03), male (1.54, 1.16–2.04), neither obese nor severely obese (1.82, 1.06–3.13 and 4.19, 2.14–8.19, respectively), had no pre-existing anxiety or depression (1.56, 1.09–2.22) or cardiovascular disease (1.33, 1.00–1.79), and shorter hospital admission (1.01 per day, 1.00–1.02). Similar associations were found in those recovering from non-COVID-19 dyspnoea, excluding age (and length of hospital admission). Interpretation Factors associated with dyspnoea recovery at 1-year post-discharge among patients hospitalised with COVID-19 were similar to those among community controls without COVID-19. Funding PHOSP-COVID is supported by a grant from the MRC-UK Research and Innovation and the Department of Health and Social Care through the National Institute for Health Research (NIHR) rapid response panel to tackle COVID-19. The views expressed in the publication are those of the author(s) and not necessarily those of the National Health Service (NHS), the NIHR or the Department of Health and Social Care. COVIDENCE UK is supported by the UK Research and Innovation, the National Institute for Health Research, and Barts Charity. The views expressed are those of the authors and not necessarily those of the funders